Ab initio shell-model calculation for ^^O in a restricted no-core model space 
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We perform an ab initio shell- model calculation for ^**0 in a restricted no-core model space, 
microscopically deriving a two-body effective interaction and introducing a minimal refinement of 
one-body energies in the spsd or spsdpf model space. Low-lying energy levels, except for the 
experimental and 2;^ states, are better described in the spsdpf space than in the spsd space. 
The structure of low-lying energy levels is discussed with an emphasis on many-particle many-hole 
states beyond the four-particle two-hole configuration. 

PACS numbers: 21.30.Fe, 21.60.Cs, 21.60.De, 27.20.-|-n 



The nuclear shell model is a successful method to de- 
scribe the structure of the nucleus [l], Q . Recently, the 
understanding of the nuclear shell structure has been in- 
creased by elucidating complicated properties, especially 
the tensor force, of underlying nuclear interaction Q . On 
the other hand, as the numerical calculation becomes 
larger and more microscopic, the meaning of the word 
"shell" has become vague. One of the ab initio meth- 
ods to describe nuclei starting with the bare nucleon- 
nucleon (NN) force in free space is the no-core shell model 
(NCSM) In the sense of using harmonic-oscillator 

(h.o.) Slater-determinant basis states, the NCSM may 
be regarded as a shell model. However, its basis states 
span a huge model space, which is taken as large as pos- 
sible, until the result of the diagonalization converges. 
Therefore, the meaning of the shell in the NCSM is some- 
what different from the usual one in the conventional 
shell model. 

While the formalism of the NCSM for heavier nuclei 
has been developed, its full application has been lim- 
ited to nuclei around the mass number A = 12 (p-shell 
nuclei), due to an explosive expansion of the Hamilto- 
nian matrix to be diagonalized. If one explores nuclei 
heavier than A — 12, some approximations or trunca- 
tions must be employed and the results obtained are no 
longer exact. However, those results may still contain 
some important physics. It is well-known that, even in a 
much more restricted model space, such as the sd or pf 
space, the conventional shell-model study with an opti- 
mal effective interaction can explain many experimental 
data 0, 0, H, Q • Thus, it is our purpose in this study to 
bridge the gap between the two approaches. 

Recently, some attempts have been proposed along 
this line. One is a shell-model calculation in the no- 
core spsdpf space combined with a minimal refinement 
of the one-body energies [l^l ■ This no-core type of shell 
model has been applied to the neutron-rich carbon iso- 
topes around ^^C, and it has been shown that low- lying 
energy levels are well described. Another is an ab initio 
shell model with a core in which a p-shell effective Hamil- 
tonian for light nuclei is constructed utilizing a NCSM 



result in a sufficiently large model space In both 

methods, no experimental information about the energy 
levels is used, and a double unitary transformation is 
performed to microscopically derive the effective interac- 
tion. In this Letter, we employ the former method and 
investigate the structure of ^^O. 

Since ^^O is a typical sc?-shell nucleus, which has the 
structure of two neutrons and a core nucleus, there have 
been many theoretical studies. As for the shell model 
with a microscopic effective interaction, the work of Kuo 
and Brown is the pioneering study [l3| ■ They constructed 
the effective interaction for the sd model space and 
showed that low-lying energy levels are well described. 
However, even in this early stage, the importance of in- 
cluding a larger space for the intermediate states in de- 
termining the effective interaction was pointed out by 
others 

In Ref. 15|, a recent sd shell-model calculation for ^^O 



using a modern NN interaction with the folded diagram 
theory has been reported. Although this method seems 
to be successful in describing low-lying energy levels, in- 
truder states having many-particle many-hole dominant 
configurations may not be properly reproduced due to 
the limited sd model space. If one wants to describe 
such many-particle many-hole states explicitly, several 
major shells must be included in the model space. In 
this sense, our present work in a restricted no-core model 
space would be preferable for the simultaneous descrip- 
tion of the ground and complicated excited states. 

Here, we outline the method of deriving the effective in- 
teraction in our restricted no-core space. The details are 
given m Refs. [13, 0. First we derive the two-body effec- 
tive interaction from a realistic NN interaction in a large 
model space Pi , solving the decoupling equation [l7| be- 
tween the model space Pi and its compliment Qi. The 
Pi space consisting of the two-body states is defined by 
a boundary number pi < 2ni + li + 2n2 + h with the h.o. 
quantum numbers {ni,li} and {ri2,^2} of the two-body 
states. Note that the full space of the two-body states 
Pi + Qi, which is necessary to accurately take into ac- 
count the strong repulsion effect of the NN force at short 
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distance reaches about 300 hft in h.o. energy. There- 
fore, in general, the diagonaUzation of the original many- 
nucleon Hamiltonian in such an energy range could not 
be performed in the foreseeable future. The effective in- 
teraction determined in this first step corresponds to a 
well-behaved interaction that renormalizes the effect of 
the short-range repulsion. 

If one uses this effective interaction for a sufficiently 
large model space in the diagonalization of the Hamil- 
tonian, this leads to the NCSM. However, this type of 
approach cannot be applied to heavier nuclei, since the 
dimension of the matrix becomes huge. Therefore, in or- 
der to construct an effective interaction suitable for a no- 
core type of shell-model calculation in a tractable model 
space, we again solve the decoupling equation between a 
smaller model space P2, such as the spsd or spsdpf space, 
and its compliment Q2 using the two-body effective inter- 
action determined in the first step, where Pi — P2 + Q2- 
Note that the boundary number pi is taken to be pi = 6 
in the present work and, thus, is not so big, although 
this number should be taken as large as possible in the 
usual sense of the NCSM. The smaller model space P2 
and its compliment Q2 are not separated by energy but 
by state, and, thus, some two-body states in the model 
space are degenerate with those in its complement. This 
degeneracy often causes a difficulty in determining the ef- 
fective interaction. Furthermore, in contrast to the usual 
NCSM, we solve the decoupling equation by introducing 
a self-consistent one-body potential for ^^O and use it 
to construct a suitable effective interaction in the small 
model space for nuclei around '^^O. This self-consistent 
procedure also makes the problem more difficult to solve, 
if the degeneracy is big. Thus, in this second step of 
the calculation, the self-consistent potential is introduced 
only for the small model space -Pj- The self-consistent 
potential is also calculated for the large model space Pi 
in the first-step decoupling. Although we decided to take 
the small number pi = 6 to make the degeneracy as small 
as possible, this artificial truncation is compensated to 
some extent by introducing a minimal refinement of one- 
body energies as explained later. 

Consequently, our shell-model Hamiltonian is com- 
posed of the sum of the one-body kinetic energy (the 
minimal refinement of one-body energies) , the two-body 
effective interaction for the P2 space, and the center-of- 
mass (cm.) Hamiltonian which is multiplied by a large 
value /3c. m. to separate the spurious cm. motion from 
the physical states [l8l |. 

In this work, we use the CD-Bonn potential [19|] as the 
original NN interaction. In Fig. [1] we present the de- 
pendence of the ground-state energy of ^^O on the h.o. 
frequency hfl. The effect of the minimal refinement has 
not yet been included in the results. The results for three 
different model spaces used in the diagonalization are 
shown: up to the four-particle two-hole (4p2h) configu- 
ration in the spsdpf space, the 4p2h spsd, and the 6p4h 
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FIG. 1: The h^l dependence of the calculated ground-state 
energies of ^*0. 



spsd from the unperturbed ground state of ^^O. The re- 
sult at the energy minimum point for the spsdpf space 
shows a less attractive energy by about 20 MeV than the 
experimental value, while the results for the spsd space 
are much more attractive. It is known that the CD-Bonn 
potential underbinds nuclei in high-precision calculations 
for few-nuclcon systems. Our result for the spsdpf space 
is consistent with this tendency. As expected, the spsdpf 
result is more favorable than the spsd results. In the 
spsd calculations, the overbinding of the energy could 
be compensated by introducing a repulsive and sizable 
three-body effective interaction, but such a calculation 
cannot be easily done. For all three types of calculation, 
there appear energy minima at around hfl = 15 MeV. 
Therefore, we choose this value as the optimal one for all 
the results in this study. 

The ground-state energy for the 4p2h spsdpf case at 
the optimal value of hfl is about —120 MeV. We note, 
however, that this number has an ambiguity due to sev- 
eral effects. In addition to some approximations and 
truncations in deriving the effective interaction and per- 
forming the diagonalization, there is a problem regard- 
ing the removal of the spurious cm. motion. Unlike 
the usual NCSM approach, our shell-model space is not 
given by energy but by state. In such a state separation, 
the influence of the spurious cm. states on the energy 
levels of intrinsic states cannot be completely removed 
by a suitable value of /3c. m.- This problem has recently 
been examined in NCSM and coupled-cluster calculations 
by Roth et ai, and they have shown that the cm. ef- 
fect is non-negligible in restricted and incomplete energy 
spaces ^2^. However, if we look at the relative structure 
from the ground state, this effect is expected to be small 
for a sufficiently large value of /3c. m., as in the study of 
carbon isotopes fio'l . 

In the calculation of the energy levels of ^^O, it is im- 
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TABLE L The excitation energies of low-lying positive-parity 
states in ^*0. The calculated results are for the cases (a) 4p2h 
spsd, (b) 6p4h spsd, and (c) 4p2h spsdpf. All energies are in 



FIG. 2: Low-lying energy levels of the positive-parity states 
in i«0. 



portant to investigate the neutron single-particle ener- 
gies for the lowest three states in ^^O. If the calculated 
single-particle levels in ^^O are not comparable to the 
experimental values, the calculation for would be 
questionable. We have found that, in the present shell- 
model space, the calculated single-particle levels are not 
so satisfactory. On the other hand, we know that, if we 
compute the single-particle energies using the CD-Bonn 
potential within the framework of the unitary-model- 
operator approach (UMOA), the results are close to the 
experiments [l^. The UMOA can be regarded as a sort 
of e'^ method or coupled-cluster theory of the Hermitian 
type 2l|. In the UMOA, the calculation of the single- 



particle energy is feasible in a sufficiently large model 
space, although it is hard to obtain complicated states. 
Therefore, in order to include the effect of the correlation 
energy of the large space into the present restricted one 
in a simple way, we change the one-body matrix elements 
of the neutron in the sd shell, utilizing single-particle in- 
formation obtained by the UMOA. This variation is the 
"minimal refinement" . 

The magnitude of the minimal refinement is deter- 
mined so as to reproduce the binding energies of the 
lowest single-particle 5/2+, 1/2+ , and 3/2+ states in 
^^O obtained within the UMOA, as a consequence of the 
shell-model diagonalization in the present space. Denot- 
ing the one-body energy for a single-particle orbit j of 
the neutron as e"'(j), we need these variations for the 
4p2h spsdpf case to be Ae"(0d5/2) = 0.99, Ae"(lsi/2) = 
-0.42, and Ae"(0d3/2) = -1.04 McV. These are added 
to the original e"(_j)'s, namely, the h.o. kinetic energies. 
We have checked that the results for the excited state en- 
ergies above the ground state are hardly changed, even 
if we introduce similar refinements into the proton orbits 
of the sd shell and the hole orbits for the neutron and 
the proton. 

In Fig. [21 we illustrate calculated low-lying energy 
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spectra of the positive-parity states for the 4p2h spsdpf 
case, together with the experimental levels. The energy 
levels in the first and second columns, labeled by "Origi- 
nal" and "Dressed" , denote the results without and with 
the minimal refinement, respectively. For the case "Orig- 
inal" , the excited levels above Ex — 3 MeV lie higher 
than the corresponding experimental ones, while those 
below Ex = 3 MeV lower than the experiments. By 
introducing the minimal refinement, the results on the 
whole become closer to experiment. A similar tendency 
was seen in the calculation for the neutron-rich carbon 
isotopes [lo| . 

To demonstrate the advantage of the calculation in the 
spsdpf space over the spsd, as in Fig. [3 the results of 
the low-lying energy levels for the three spaces in the 
"Dressed" case are given in Table HI Note that the mag- 
nitudes of the minimal refinement of the one-body en- 
ergies depend on the model-space size employed. The 
magnitudes for the 4p2h spsdpf have been given be- 
fore; those for the 4p2h spsd are Ae" (0^5/2) = 2.12, 
Ae"(lsi/2) = -3.14, and Ae"(0d3/2) = -1.81 MeV, 
and for the 6p4h spsd, Ae"(0d5/2) = 2.42, Ae"(lsi/2) = 
-2.59, and Ae"(0d3/2) = -1.30 MeV. These values are 
determined in each restricted model space so as to re- 
produce the UMOA results. The energy spacings for the 
4p2h spsd case are rather compressed compared to the 
experimental data. The inclusion of the GpAh configura- 
tion does not cause any significant change in the relative 
spacings. The results for the 0^ state for the spsd cases 
are fairly close to the experimental rather than the 



experimental Oj. However, by including the pf shell, 
the calculated state becomes close to the experimen- 
tal 0;^, and the energy spacings are expanded. Then the 
results become, on the whole, close to the experimental 
values. This suggests the risk of using the experimental 
single-particle energies in the conventional shell model, 
leading to an accidental reproduction of energy levels of 
the wrong structure [2^ . 

As shown before, the minimal refinements for the spsd 
cases are considerably larger than the spsdpf ones. This 
means that the minimal refinements in the spsd cases 
are not so realistic. In general, the effect of the many- 
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TABLE II: The excitation energies in MeV and the prob- 
abilities of the 4p2h configuration P4p2h for the lowest 4p2h- 
dominated 0^ state in ^**0, for the configurations up to 4p2h, 
6p4h, and 8p6h in the spsd space. 



Configuration (spsd) 
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8p6h 
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body effective interaction becomes bigger as the size of 
the model space becomes smaller. Therefore, in the spsd 
space, the change of only the one-body energy may be 
too simple a treatment. In such a small model space, we 
should change not only the one-body energy but also the 
two-or- more-body matrix element. 

Interestingly, in the spsdpf case, the result of the 
state appears at near the experimental 0;^, and there is 
no calculated 0"*" state near the experimental Oj. Fur- 
thermore, we cannot find a 2"*" state near the experimen- 
tal 2^. These experimental and 2J^ levels have been 
treated as (deformed) 4p2h states in shell-model calcu- 
lations 1221, while these states have been considered to 



be well-developed ^'^C+a cluster states in Refs. 23l. l24|. 
In our calculation, the 4p2h configurations are included. 
Nevertheless, there occurs no indication that our 4p2h 
states come down to the low-lying energy region. 

There may be some reasons for the disappearance. One 
is the simple treatment of the minimal refinement, even in 
the spsdpf space. Another is the effect of higher config- 
urations, such as the 6p4h, on the 4p2h dominant state. 
To investigate this effect, we have calculated the lowest 
4p2/i-dominated state in the spsd case. In Table HIl we 
tabulate the excitation energies and the probabilities of 
the 4p2h configuration for the lowest 0"*" state dominated 
by the 4p2h configuration, for the configurations up to 
4p2h, 6p4h, and 8p6h in the spsd space. The probabil- 
ities of the 4p2h configuration contain the contributions 
from all the 4p2h states in each space. The calculated 
4p2h state is the fourth 0+ state in the spsd case. The 
third 0+, which is dominantly composed of the (0^3/2)^ 
configuration of the neutron, lies at about 16 MeV in ex- 
citation energy for the three cases. Although the excita- 
tion energies in TableHIlai'e far from the experimental 0^, 
a considerable contribution of about —10 MeV from the 
6p4h configuration is brought to the excitation energy, 
while that from the 8p6h configuration is rather small. 
This indicates the importance of also including the 6p4h 
in larger spaces, even if the 4p2/i is the dominant config- 
uration, when (spherical) h.o. wave functions are used 
as the basis states. Such a calculation is a challenging 
problem to reveal the real structure of the experimental 
O2 and 23 states. 

In conclusion, low- lying states in ^^O are well described 
by our ab initio shell-model calculation with minimal re- 
finement in the no-core spsdpf space. However, a larger 
model space is necessary to properly describe the exper- 



imental 0^ and 2j^ states. 
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